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SUMMARY 


Five quasars are chosen which all show evidence for the presence of intervening 
galaxies. An analysis of the observational data results in one of these— 
3C 268.4— being a good candidate for the occurrence of a visible secondary 
gravitational lens image. The expected apparent (V) magnitude of this 
companion image is 21:2* 19, and if confirmed by observations suggested, 
would provide an independent measure of the Hubble constant. 


I. INTRODUCTION 


The idea that a ray of light might be bent by the gravitational action of a massive 
body was first suggested by Newton in his Opticks (1704). However, it was not until 
more than two centuries later in 1919 that the effect was observed. In that same 
year a related phenomenon, namely the gravitational lens effect, was first described 
by Sir Oliver Lodge—an effect which to date has still defied observation. In 
addition to bending the light from a background source, a fainter secondary image 
of the source is produced on the opposite side of, and much closer to the gravita- 
tional deflector. As the alignment of the source, lens and observer is improved the 
two images of the source merge into a ring image, resulting in an extremely high 
amplification in apparent brightness of the background source. This high degree 
of alignment is however most unlikely, and in this paper the observational evidence 
for the production of a weak secondary lens image of a background source by a 
gravitational deflector is investigated. 

Since their discovery in 1963 quasars have been associated with the gravita- 
tional lens effect both in the role of background sources and as amplified images of 
much fainter objects (Barnothy & Barnothy 1968; Press & Gunn 1973; Jaroszynski 
& Paczynski 1974; Gott & Gunn 1974). As the former possibility is far more likely, 
five quasars—3C 196, 286, 287, 268.4, PHL 938—which for various reasons exhibit 
evidence of intervening galaxies, were chosen as possible candidates for lensed 
background sources. 


2. EVIDENCE FOR INTERVENING GALAXIES 


The most direct evidence for the presence of an intervening galaxy is if its 
redshift, zp, is low enough (2p 0:5) to render the galaxy visible photographically 
(magnitude <23™). In addition the galaxy would have to be separated from the 
quasar by at least ~ 2 arcsec, otherwise, the galaxy would be thought of either as 
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an underlying galaxy in which the quasar was embedded, or a jet associated with 
the quasar. Indirect evidence—where no galaxy is seen photographically—tends 
to rely on the effects of the gas associated with galaxies and hence implies the 
presence of a spiral, rather than an elliptical intervening galaxy, since the latter 
contain very little gas. In general, spiral galaxies are not compact enough to act 
as gravitational lenses (Sanitt 1971; Jaroszynski & Paczynski 1974), and a difference 
in optical absorption of up to ~2™ is expected between the two lensed images 
(Wagoner 1967). Therefore it is unlikely, though not impossible, that a visible 
secondary lens image is present in cases of indirect evidence for intervening 
galaxies. | mE 

Three indirect effects considered are 21-cm absorption, Faraday rotation and 
optical absorption lines (Wagoner 1967). All these are evidence of the presence of 
an intervening galaxy, although the latter two effects could also. certainly be 
intrinsic to quasars. It should also be borne in mind that since the lens effect is 
independent of the wavelength of e.m. radiation, compact radio components of 
quasars may also be lensed by intervening galaxies, producing secondary radio 
lens images. Before examining the observational evidence for each quasar separately, 
the various equations and assumptions which enable estimates to be made of the 
likelihood of a secondary lens image being visible are considered. 


3. ASSUMPTIONS 


We consider an axially symmetric elliptical galaxy acting as a gravitational lens 


whose mass distribution follows the de Vaucouleurs (1948) law, with a constant 


mass : light ratio throughout the galaxy. 
The projected surface density p is given by 


logio 2 = —3:33(a1/4— 1) (1) 
pe 


where « = d/re, d is the distance to the centre of the galaxy, and the subscript ' e’ 
denotes quantities evaluated at the ‘ effective radius ' (re) of the galaxy. (This is the 
radius which contains half the total mass of the galaxy.) The mass M(d) within a 
cylinder of radius d and axis along the line of sight to the observer is given by 


M(d) = f 2mxp dx. - : | (2) 


For simplicity we assume (a) a constant mass : light ratio, and (b) axial rather 
than elliptical symmetry. Assumption (a) is justified since no definite statement can 
be made, through lack of data, on how the mass : light ratio of elliptical galaxies 
varies with radius. Also, since we are mainly concerned with the mass distribution 
from ~0-5 to 20 kpc, the existence of a massive halo would hardly affect the 
results. | 

The lens effect for elliptically-symmetric deflectors has been considered by 
Bourassa, Kantowski & Norton (1973). However, in the present context, no gain in 
accuracy of the predictions in this paper can be obtained by such a more high- 
powered mathematical model, which takes into account elliptical symmetry, in 
view of all the inherent uncertainties of the mass model for galaxies employed. 
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4. CALCULATIONS 


The apparent positions of the primary and secondary images with respect to the 
centre of the lens galaxy (d) are the solutions (Sanitt 1971) of 


_ 4G (2s— 29) (1-450) a (90) pyi p Ben (00) d 
diii wv Ten, M(a)+ gu — (3) 


where G, Ho, c are the gravitational constant, Hubble constant (adopted throughout 
as 50 km s-t Mpc~), and velocity of light; zy and zs are the redshifts of the galaxy 
and quasar respectively; B is the alignment angle between the galaxy and unlensed 
source position; cz'(25)/ Ho is the ratio of the intrinsic to angular size of a rigid 
rod at zy, 


a (2v) = un re [(x 4-29) — (1 - 2)! 7] (4) 


(qo = 4, A = o Friedmann model); T (2p, 2g) is a cosmological correction factor 
which is close to unity (tabulated by Refsdal 1966). 

Equation (3) can be solved by determining the intercept points of the following 
two curves 


Ya = f(a) (5) 


2 (Lue j| o MocanzsT (m zs) [P a | Bez (2p) | 
Y» = (reus [iG an) eM fate E (6) 


where F signs refer to the primary and secondary images, and f(x) = M (a)/Mtot- 
From (1) and (2) 


f a’ exp ( — a(o/1/4— 1) do’ 
— Jo 


JC. cause a (7) 
f a’ exp {—a(a’l/4—1)} da’ 
where a = 7-669 (o is a dummy variable). 
Therefore 
fla) = 1 [expan] {X> S atl 8) 
or 
fla) = [exp- as] | È it (9) 


The equations (5) and (6) do not necessarily always have a solution. If we 


define Bo as 


Mio ([4G(2s — zp)(1 + 2p) da) 
Po = 0°73 f* | | CHzyzs T (2, Zs) (10) 
then provided 
p EY (typically = o-1 arcsec) | (11) 


a secondary image will not be produced if 


B> Bo (i.e. no solution to (5) and (6)). 
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Po is defined in terms of the de Vaucouleurs law. However, other mass distributions 
exhibit similar behaviour, in that a primary image situated more than a given 
separation from the galaxy, does not produce a secondary image. It was effectively 
this constraint that Gott & Gunn (1974) applied to the case of the double quasar 
1548 4- 115, which turned out not to have a secondary lens image down to a limiting 
magnitude of 23?! (Gott & Gunn, private communication). 
The relative intensity of the primary and secondary images, R, which in the 
case of a point mass is oj/o2, is modified, in the case of a distributed mass, by a 
‘thickening factor’ which arises from the relative gradients of the two functions 


(5) and (6) (Sanitt 1974). Thus 


 dY4 [dYs 
R = ai dx [| da |a : (12) 
as de dYA d Yg 
da | dx |a 


Equation (6) contains the factor re/Mo91!/?. A mean value and standard deviation 
of 


re kpc 
Mair = F9tio7 | (13) 


is obtained from a list of 25 elliptical galaxies (Fish 1964), and adjusted for a 
Hubble constant of 50 km s-1 Mpc7! (Mi is Mio in units of 1011 Mọ). 

In general the angular separation between the secondary image and intervening 
galaxy will be less than 1 arcsec, so that the two images will be unresolved. In 
order to separate out the individual contributions of the galaxy and secondary 
image to this combined image, it is necessary to relate the mass of the galaxy—which 
determines the apparent magnitude of the secondary image Msec, to the apparent 
magnitude of the galaxy mga). The two magnitudes msec and mga together give the 
apparent magnitude of the combined image, "comp, which is observed. We use 
the relationship between the mass and (luminosity)?/? (Noonan 1971) 


Mi | 
+0-101 
1o 06M 110 ~ "075 —0-052° (14) 


This corresponds roughly to a mean mass : light ratio for elliptical galaxies of 
^19 (Ho = so km s! Mpc“!). The apparent magnitude of the galaxy, Mgal, is 
then determined using K-corrections K (zy) from Schild & Oke (1971), 


Mgal = My + K (ap) + 43:88 + 5 logio 2[(1 + 2) — (1 + 25)!/2] (15) 


(a qo = 4, A = o, Friedmann model is used). 

The aim of the calculation is to determine the apparent magnitude of the 
secondary image—Msec, and the galaxy—#mgal, assuming that the two images are 
unresolved and only the combined magnitude mcomp is known. 

In addition other measured parameters are 2p, Zs and 0"—the separation between 
the primary image and the mean position of the combined image. The following 
iterative procedure is used. | 

(1) The difference in magnitude between the primary (prim) and secondary 
images, and the quantity (B/re) are both considered free parameters. The dimension- 
less distances from the galaxy to the primary and secondary images a1, o2 can then 
be determined. 
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(2) A trial value of the apparent magnitude of the galaxy—ga) is assumed, 
and the mean (dimensionless) position of the combined image (weighted by 
apparent magnitude) can be determined. 

(3) The scaling factor re can then be determined since 6” is known. 

(4) The relationships (13), (14), (15) enable the magnitude of the galaxy— 
Mga) to be computed from re. 

(5) This value of mga, is compared with the trial value in (2) and the process 
from (2) to (5) is repeated until a consistent value of mga is determined. By varying 
the parameter (B/re), msec as a function of Mcomp can be generated together with all 
the other unknown parameters. (Note that the relationship between msec and 
Mcomp and other dimensionless parameters are all independent of Ho.) 


5. INDIVIDUAL CASES. 
3C 196 
This quasar has the highest redshift (0-871) of a group of four quasars which 
Kristian (1973) lists as having close off-centre components. The seeing limited 
nebulosity (~23™) 1-5-2 arcsec SE of the quasar is more likely to be a jet rather 
than an underlying or intervening galaxy. However, even if the latter possibility 
is the case, a secondary image would be below 23™ and thus too faint to be visible. 


3C 286 


No companion image to this quasar is visible on a 200-in. photograph (Kristian 
1973) (Mcomp 223™). However, a 21-cm hydrogen absorption line at a redshift 
Zb = 0:692 (Zem = 0:849), is present for which the most likely explanation is an 
intervening galaxy. Application of the calculations in Section 4 together with 
Msec %23™ (Mprim = 17:3), yields an upper limit of — 101? Mo for the proposed 
intervening galaxy. The radio source is single (diameter ~ 0:03 arcsec, Broten et al. 
1967). Since no secondary radio image is visible with the 5-km telescope (Pooley, 
private communication), then the following limits can be obtained: 


(1) If the secondary image is resolved from the primary (separation 2 3 arcsec); 
then the detectability limit (12-hr run, good tropospheric conditions) is ~o-o1 Jy; 
1.€. Msec Z 2411-3. 

(2) If the secondary image is unresolved, then the detectability limit is 0-12 Jy 
(~2 per cent of the primary image flux) obtained from limits on the variation of 
visibility with hour angle. Since this implies msee 221™-6, the optical limit of 
Msgec 2 23™ is more important. — 


Thus even down to 24™ it is quite probable that no secondary image would be 
visible. 


3C 287 


This quasar (my = 17-67) has a companion 3:4 arcsec away and ~ 3™-5 fainter 
(blue) (Matthews 1968; Kristian & Peach 1968). There may also be some faint 
nebulosity although this is very doubtful. The radio source is single and has a high 
rotation measure (Donaldson, Miley & Palmer 1971). In order to fit the data to 
that of a combined galaxy-secondary image, the intervening galaxy must have 
Zp 20°5, Miot ~ 2x 101? Mo and msec ~ 21:2. Since no secondary radio image 
is visible on the Cambridge 5-km telescope (Pooley, private communication), this 
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implies mgec Z; 2271-239? (intensity of source = 2'4 Jy); since the two values of 
Msgec do not agree, the assumption that the observed companion is a combined 
galaxy-secondary i image is false and therefore the xcd is not producing a secon- 
dary lens image. 


PHL 938 


This quasar represents the most likely case of an intervening galaxy causing 
absorption lines (Zapg = 0:613, Zem = 1:955). Applying the calculations described 
in Section 4 and assuming a galaxy mass of say 1011 Mọ and msec = 23™, we obtain 
a primary separation of ~2 arcsec (1477 kpc) (qo = 4, A = o model). Since the 
separation could be as high as ~20 kpc and still result in quasar absorption lines, 
it is not clear whether this quasar would show a visible secondary image or not; 
however, a photograph of this quasar down to a limiting magnitude of ~23™ 
might exhibit a companion image. 


3C 268.4 


A 200-in. plate shows a faint (~20™-5) companion image 2:5 arcsec away in 
position angle ~ 170? from this quasar (mprim = 18-42) (Kristian 1973). Kristian 
also reports a high density of galaxies which indicate a cluster at zy ~ 0°35 plus 
some foreground galaxies. Assuming that the companion is a combined secondary- 
galaxy image, the latter having a redshift —0:35, the calculations of Section 4 
result in msec = 21:2, i.e. the secondary image accounting for just over half of the 
brightness of the combined image. The absence of optical absorption lines would 
tend to indicate an elliptical galaxy (no gas), and the calculations result in a mass 
of — 5:7 x 101! Mo. Fig. 1 shows the full results of the calculations; also indicated 
are loci of constant time delay between the primary and secondary images (see 
Appendix). 

The derived parameters are (see Fig. 2): 

Intervening galaxy (Ho = 50 km $1 Mpc) 
Miot = 5:65 x 1011 Mo 
My = —21'5 
Mgal = 21°3 
7e = 074 arcsec = 4:45 kpc 
(Mass/light)y = 17:8 
Bo = 2:2 arcsec 
maximum observed separation of primary which produces a secondary image = 
2:9 arcsec. 


Lens images 
B = 1:46 arcsec 
Msgec = 21°2 


dprim = 2:26 arcsec = 13:6 kpc 


dgec = 0°45 arcsec = 2:7 kpc 
Atg = 110 days 

Atp = 165 days 

Attot = 275 days. 
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Fic. 2. Positions of radio components + , and optical images of the primary—P, secondary— 
s, and galaxy—G. Unlensed positions are circled. 


Uncertainty in msec 
The contribution to the uncertainty in msec for each of the assumptions and 
empirical relations are as follows: 


Contribution to 
uncertainty in mygec 


Assumption (m) 
(1) Cosmological model—a range of + o'5 in the value +Or1 
of go | 
(2) Apparent magnitude of meomp + 05 +1°0 
— [e] ° 9 
(3) Error of + 0-05 in zy (assuming that zp ~ 0°35 is | —0'3 
appropriate to the value of the distant cluster) | 10:5 
(4) Uncertainty in mprim due to variability +o™-5 - to:06 
(5) If a point mass rather than a de Vaucouleurs model —0'3 
is used | 
(6) re Mto! = 11907 +0°2 
| —O*I 
(7) Ratio Mi1/1070 6My-11-0 = 0-075 --o- ror —0'4 
—0*052 T 0:75 
Combined uncertainty 1:6 
—1'2 
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Observed variability in 3C 268.4 (Peach 1969) would be mimicked in the 
secondary image and delayed by ~ yr (Ho = 50 km s^! Mpc71). A measure of 
this time delay would result in an uncertainty of ~ +28 per cent in Ho, dropping 
to ~ +14 per cent if 7comp were accurately T "The time delay is reduced 
6:2 per cent if the qo = 1, A = o Friedmann model is used (Atg = 105 days, 
Atp = 153 days), and increased by 5:4 per cent for the go = 3, ‘ clumpy ' model 
used by Gott & Gunn (1974) (Atg = 126 days, Aty = 164 days). 


Radio structure 


A radio map of 3C 268.4 (Pooley & Henbest 1974) shows an asymmetric 
triple source, with the weakest middle component (0:05 + 002 Jy) coincident with 
the optical image (see Fig. 2). The position of this image (128 06™ 428-11 + 08-03, 
43^ 56 2":1 € o":3) is a weighted mean position of an optical position measured 
by this author and corrected for the effect of the faint companion image, and the 
radio position quoted by Pooley & Henbest (1974). A secondary radio lensed 
image would have a flux —0o:004* 9005; Jy, which is below the detectability limit 
of the 5-km telescope (0-o1 Jy for a 12-hr run with good tropospheric conditions). 
Four days observations, however, would lead to a factor 2 improvement in signal/ 
noise, which could be sufficient to yield a positive measurement. The displacement 
of the optical and middle radio component positions would be — 0:8 arcsec, and 
0:5 arcsec for the strong radio component (0°57 Jy) which is 3:9 arcsec from the 
galaxy. This radio component is too far away from the galaxy (> 2:9 arcsec) to 
produce a secondary image. 


6. DISCUSSION 


An important question is what alternative explanations are there for the com- 
panion image to 3C 268.4, other than what has been discussed so far? 'There are 
basically three possibilities: 


(1) It 1s a jet associated with the quasar. 

(2) It is a galaxy but not associated with the cluster at x = 0°35. 

(3) It is a member of the cluster at z = 0:35, but no secondary image is visible 
down to 23™. (This is basically the case for the double quasar 1548 + 115 considered 
by Gott & Gunn 1974.) 


Since a 200-in. photograph taken by Kristian (private communication) shows 
that the companion image is resolved from the primary and does not appear 
connected, (1) is an unlikely possibility. 

The proximity of two nearby galaxies is such that (2) is alse unlikely. 

If (3) is the case then the mass : light ratio of the galaxy must be much lower 
than normal. A value of —5-7 would be indicated compared to a mean value 
assumed of ^ 19. Out of 13 non-peculiar elliptical and SO galaxies considered by 
Noonan (1971), only two have a mass : light ratio as low as 5-7. 

Thus if (3) is operating the galaxy must have an abnormally low mass : light 
ratio, and also an upper limit to its mass of ~4 x 101! Mo. 


7. CONCLUSIONS 


Out of five quasars considered which exhibit evidence of intervening galaxies, 
one case—3C 268.4—should show, subject to reasonable assumptions, a visible 
secondary lens image. 
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Observational proof would be possible if a spectrum of the combined image 
showed emission lines similar to those in the main quasar, if variations in brightness 
of the main quasar were duplicated by the combined image, and if the companion 
image were shown to be double. 

The spectrum of the combined image would consist of an approximately equal 
mixture of a galaxy at a redshift of ~ 0:35, and a quasar at a redshift = 1-4. The 
strong quasar emission lines of C 1v and C 111] would be visible at 3722 and 4578 A 
(Schmidt 1968), and the galactic H and K lines would be at 5357 and 5311 A 
(2p = 0°35). 

3C 268.4 has been observed to vary by ~o™-5 in 1 yr (Peach 1969). The 
combined image would therefore also vary ( ^» 01-25 variation) with a time delay of 
~275 days (Ho = 50 kms! Mpc~'). However, because of the proximity of the 
combined image to the primary, variations in brightness of the combined image 
will be difficult to measure. 

The separation between the galaxy and secondary image would be 0:45 arcsec; 
however, the images could only be resolved under the most superlative of seeing 
conditions of the 200-1n. telescope. 

Finally, it is suggested that a photograph of PHL 938 taken down to a limiting 
magnitude of ~23™ might exhibit a visible secondary lens image. 
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APPENDIX 
TIME DELAYS 


The two lensed images of a variable source will exhibit identical variations in 
brightness, but with a relative time delay due to differences in (a) path length, and 
(b) gravitational potential along the two rays of light. Following Cooke & Kantowski 
(1975) the geometrical time delay due to (a) 1s 


D | 
cAtg = $(1+2p) = —— [(dsec + x)? — (dprim — x)?] (1) 
DysDp 


where D», Ds, Dvs are the ‘ apparent size ' distances from the observer to the galaxy 
and source, and from the galaxy to the source, which can be found in Refsdal 
(1966); and x is the projected distance from the centre of the galaxy to the position 
of the unlensed source. The potential time delay Atp is 


_ 4G dorim M (d) | 
ety = ST (remp) | PE) aa (2) 
i.e. 
cAtp days = 22:87 (1 + zy) Mir [5 do (3) 
ag Q 
substituting « = yt and with a = 7:66925 we have 
“$day (ert cay 
J. c ds 7 7 " alit y m dl dy | (4) 


ap a aj4 2 7 
= nil itd ay | e | Am dy) (5) 


alh y alst 
where we have used 
7. gi 
E 4 ^ 
f(x) = 1x — exp — aa [X j^ | (6) 
see equation (8) main paper. 
The integral 


ai e—ay 
dy 
an y 


can be expressed in terms of the exponential integral 
oo e` 
E(x) = I u, 

zs x 


and the second integral is 


ec y y eX (7) 
7 i=0 j=i+1 11j 
us 
al . 
[7/9 da = 1n 91— 4 Esos!) - Erao) + Afo) A t] (8) 
ag 
where 


§ 
Aene E (9) 
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Fic. 3. Contours of constant I (a1, «2) in the «1, «2 plane, see text. 


with 
bo = 2:59286 ba = 73:4452 
bi = 1221602 bs = 68:43483 
b2 = 32°13950 bg = 40:37260 
bs = 57:10166 
For large x, Ei(x) ~ 1/x e? or more precisely for x> 1(a>0-0003), 
x2 4- 3x 4- Lo 
z e LUN 
x e* Ex) "Oc (10) 
where | 
li = 2334733 mı = 3330657. 
ly = o:250621 mz = 1°681534. 


The Fj(x) terms have only been included for completeness, and are in fact 
negligible in normal ‘ lens ' situations. 


The In o;/o2 term is that due to a point mass (of mass Miot) and it is con- 
venient to write the time delay in the form 


ctp = £5 + Zp) Mio In ^5 Ion ag). (11) 


This defines the correction factor J (o; a2) which represents the departure from the 
point-mass approximation. 
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Fig. 3 shows contours of constant (o1, o2) in the a1, «2 plane. Since a ' distri- 
buted ' mass is a less efficient gravitational lens than a point mass, the time delay 
is less and Z(o3, o2) « 1. As expected as o3, ag > oo, the galaxy behaves more like a 
point mass and Z(o, o) — 1. The function can be represented reasonably accura- 
tely by the approximate formula: | 


I(o3, ag) = I—exp {— 0*7(a30)9/8). (12) 


This gives values of I(o1, o9) accurate to 10 per cent for o:15 « ojxg « o*5, and 
3 per cent for ajag>0°5. 
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